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Abstract. We calculate the magnetization profiles of the a x and crj operators for the XX- 
model with hermitian boundary terms. We study the profiles on the finite chain and in the 
continuum limit. The results are discussed in the context of conformal invariance. We also 
discuss boundary excitations and their effect on the magnetization profiles. 



1. Introduction 

We continue our study of the XX-model with boundaries defined by the Hamiltonian 

x L-l 1 

H = - 22 [o+o- i+1 + o i <t+_i] + -= [a-ff! + a+af + a z cr z + + /?_a L + (3 z a z L ] (1) 
i=i 

where af = |(a"J±icx|). By a x , a y and a z we denote the Pauli matrices. Since we restrict 
ourselves to hermitian boundary terms, the values of a z and f3 z are real numbers and a + 
and (3 + are the complex conjugates of a_ respectively In the first publication Jl|] of 
this series we have shown how to diagonalize H by introducing an auxiliary Hamiltonian 

L-l 



+— = [a-o%<r x + a+a x Q a+ + a z a z + (3 + ata x L+1 + (3_a L a x L+1 + f3 z a z L ] (2) 
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(following an idea of 0) which in turn may be diagonalized in terms of free fermions ||. 
Note that H\ ong commutes with erg and cr£ +1 . Hence the spectrum of H\ ong decomposes 
into four sectors (+, +), (+, — ), (— , — ), (— , +) corresponding to the eigenvalues ±1 of 
Uq and cx£ +1 . The spectrum and the eigenvectors of H are obtained by projecting onto 
the (+, +)-sector. 

In this paper we are going to consider the profiles of the spin operators a z and a x . 
We have already shown in ]TJ how to compute the one point function of the crj-operator 
in terms of pfaffians. The <r|-profiles are easier to compute. We will study the ground- 
state profiles as well as the profiles of certain excited states of H. The excited states of 
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H we are going to consider correspond to single fermion excitations with respect to the 
ground-state of Hi ong , which are characterized by a non- vanishing excitation energy as 
L goes to infinity and which correspond to boundary bound states. 

We are able to obtain exact results for the ground-state profiles only for certain 
values of the boundary parameters In the case of the cr|-profiles we will find exact 
expressions on the finite chain, in the continuum limit L ^> 1, z — j/L fixed and in 
the limit L — > oo, j ^> 1. The cxj-profiles are given in terms of pfaffians of 2j x 2j 
matrices, that we have not been able to compute on the finite chain or in the continuum 
limit exactly, therefore our exact results for the crj-profiles are restricted to the limit 
L — > oo, j > 1. 

A combination of extensive numerical computations (for lattice lengths L ~ 600 ) 
and of our exact results has led us to conjectures for the profiles of <r| and crj in the 
continuum limit for more general values of the boundary parameters. Note that once 
we know the <7?-profile for one specific choice of the boundary terms, we know already 
the profile of cxj for the boundary terms with a x and a y interchanged in ([!]). 

We will discuss our results for the ground-state profiles in the context of conformal 
field theory. According to Burkhardt and Xue 0, in the continuum limit z = 
j/L,L 3> 1 the general form of the profile of any scalar scaling operator with bulk- 
scaling dimension x^ is given by 

(<P( 3 )) ab = [(L/tt) sm(7rz))- x *F ab (z). (3) 

The functional form of T a b{z) depends on the boundary conditions a and h at the 
left respectively the right boundary. We have shown in || that the XX-model with 
hermitian boundaries corresponds to the free compactified boson on a cylinder with 
Neumann-Neumann, Dirichlet-Neumann or Dirichlet-Dirichlet boundary conditions. 
The bulk-scaling dimensions x^ obtained from the two-point functions for the periodic 
chain are 1 respectively | for a z - respectively o x - 0. . 

Conformal field theory makes more restrictive predictions for the behaviour of the 
profiles of scaling operators near the boundary. From (|3]) we obtain that near the 
boundary the profile is given by 

(<P(j))a b =r X ^a b (z) *<1. (4) 

It was already argued in that if the scaling function ip a b(z) is different from zero near 
the boundary, it has the asymptotic form 

ip ah (z) = A[l + Bt b z 2 + • • •] Z <1. (5) 

Furthermore it has been shown in 0, || |l0| that the amplitude B^ b is related to the 
Casimir amplitude A a b via 

Bt b /A ab = -Anxo/c (6) 

where c denotes the central charge of the theory. The Casimir amplitude A a b is defined 
by the large L behaviour of the ground-state energy £<g, i.e. 

E = eoo L + f ab + AabL- 1 + ■■■ (7) 
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where by e ro and f a b we denote the bulk respectively the surface free energy. The 
Casimir amplitudes for the boundary terms of the Hamiltonian (|1]) can be read off from 



the partition functions obtained in H, I I 



As we are going to show, when dealing with boundaries, it is not a priori clear which 
combinations of spin operators correspond to the scaling operators of the continuum 
theory. Our conjectures for the continuum limit of the profiles will yield, that the 
profiles of aj and of cr| indeed have the form of (§) (see however equation (|35l)). 

We will see that if the boundary conditions we are going to consider correspond 
to Dirichlet-Neumann of the free boson field, the profiles satisfy (d) and (P). We will 
also see that for the boundary terms corresponding to Dirichlet-Dirichlet and Neumann- 
Neumann boundary conditions the profiles show a more complicated behaviour near the 
boundary. In these cases the equations ([D and (|6]) do not hold in general for the profiles 
of (Tj and cr|. 

In the case of Dirichlet-Dirichlet boundary conditions we have to build an 
appropriate linear combination of cr| and cr| +1 in order to obtain the correct behaviour 
of a scaling operator. However, in the generic case the linear combination we are going 
to obtain works only at one boundary at a time. The cr^-profile vanishes exactly in this 
case. 

For the case of Neumann- Neumann boundary conditions we obtain the correct 
behaviour by considering a linear combination of crj and crj. Again this works only at 
one boundary at a time in the generic case. Here the aj-profile vanishes if no diagonal 
boundary terms are present at the boundary (i.e. a z = f3 z = 0). However, we have 
also studied a type of boundaries corresponding to the Neumann-Neumann boundary 
conditions, where diagonal and non-diagonal boundaries are present at a boundary at 
the same time. In this case the profile of the cxj-operator vanishes only in the leading 
order and is determined by a secondary field instead by a primary field. 

Our results in the limit L — > oo,j ^> 1 are in agreement with the results of 



Affleck [12]. He has considered the semi-infinite XXZ spin-1/2 chain with a transverse 
magnetic field applied at the end of the chain (this corresponds to a + = a_ 7^ 0, 
a z = /3 Z = (3 + = P_ = and L — > 00 in our notation). Affleck has used the bosonization 
technique, whereas our results have been obtained on the lattice. Assuming Neumann 
boundary conditions for the boson field, he obtained the cxj-profile for large values of j. 
At the free fermion point the result reads 

(of) = c-iycr* (8) 

where the constant C does not depend on the strength of the field and at the free fermion 
point its numerical value is given by 

C = 0.912 171 210446 (9) 

The case of a longitudinal magnetic field at the end of the chain has also been considered 
in fl2| (corresponding to a z 7^ 0, a + = a_ = j3 z = (3+ = = and L — > 00). Imposing 
Dirichlet boundary conditions for the boson field, the large j behaviour of the cx|-profile 
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has been evaluated for small values of a z , yielding at the free fermion point 

(of) « Ba^-iyr 1 (10) 

where the value of B is undetermined. The case of diagonal and non-diagonal boundary 
terms at the boundary at the same time has not been considered. The boundary 
conditions where we have been able to compute the profiles exactly in the limit 
L — ► oo, j ^> 1 can always be mapped onto the boundaries treated by Affleck. Our 
results obtained on the lattice agree with (§) and fllPl) . For the cr|-profiles we have not 
restricted ourselves to small values of a z and we have derived the a 2 -dependence of B. 

As already mentioned above, we will also consider the magnetization profiles of 
certain eigenstates of H, which are obtained by the excitation of a massive fermion with 
respect to the ground-state of H\ ong . The existence of such fermions has already been 
mentioned in for special values of the boundary parameters. Before considering the 
profiles, we are going to study the appearance of such fermions for general, hermitian 
boundary terms. 

We will see that massive fermions appear if the values of certain functions of the 
boundary parameters become greater than a certain threshold. The respective excitation 
masses will turn out to depend only on the boundary parameters of one boundary at a 
time. In the general case (non-diagonal together with diagonal boundary terms) we will 
find at most one mass per boundary. 

In the case of purely diagonal boundary terms the fermion masses correspond 
exactly to the non-vanishing energy gaps obtained in the Bethe Ansatz by comparing 
the energy of the reference state with the energy of the one-magnon excitations, which 
correspond to boundary 1-strings |[l3| . 

We will see, that the magnetization profiles for the excited states under study 
here, show an exponential fall-off into the bulk for diagonal as well as for non-diagonal 
boundary conditions (the profiles for the non-diagonal boundary terms have been 
obtained numerically). Hence we are going to refer to these states as to boundary 
bound states. These profiles appear also for H since the excited states we are going to 
study lie in the (+, +)-sector of the Hilbert space of H\ ong and do therefore correspond 
directly to eigenstates of H. 

This paper is organized as follows: We start with the exact results in section 2. 
These results are restricted to certain choices of boundary terms. In section 3 we will 
give our conjectures for the profiles in the continuum limit for more general types of 
boundaries. We are going to discuss these results in the context of conformal field theory 
in section 4. Boundary bound states are the subject of section 5. Our conclusions will 
be given in section 6. The details of our calculations will be given in the appendix. In 
Appendix A we shortly discuss the diagonalization of H\ ong and the projection onto the 
(+, +)-sector in order to provide the basic facts and notations being necessary to follow 
our calculations. Appendix B deals with the exact computations of the crj-profiles. How 
we arrived at our exact results for the crj-profiles is shown in Appendix C. In Appendix 
D we discuss the numerical verification of our conjectures of the ground-state profiles in 
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Table 1 . The values of the boundary parameters where we have computed the ground- 
state profiles of crj exactly on the finite chain and in the continuum limit. By DN 
(Dirichlet-Neumann) and DD (Dirichlet-Dirichlet) we denote the respective boundary 
conditions on the free boson . For case e there appears a free, real parameter 8. 
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the continuum limit. 
2. Exact results 

For technical reasons explained in |Appendix A] , we have been able to obtain analytical 
results only for certain choices of boundary terms (see table |l] and table ||). We will begin 
with the results for the crj-profiles on the finite chain, before turning to the respective 
expressions in the continuum limit. Thereafter we will give our exact results in the limit 
L — > oo, j ^> 1 for the profiles of crj and crj. 

2.1. Exact results on the finite lattice ({crfj- profiles) 

We have computed the exact profiles of aj for values of the boundary parameters given 



in table [T] [|. Note that for technical reasons (see Appendix B for details) our result for 



case e is restricted to the choice 6 = in this section. Furthermore, observe that the 
ground-state is twofold degenerate for case d with L odd and for case e with L even, 
due to the existence of additional zero modes 0. In these cases, we have obtained 
the profiles for the two ground-states, which have a well-defined fermion number (see 



Appendix B|) . These profiles are indicated by a ±-sign in the respective expressions. 
• Case a 

(<x?) = ( (-1) L tan + (-lVcot ( 2j -^ n ) . (ii) 



Case b 



Case c 



w = izi)i cot %zJ> (12) 




-r 



«) = ^l(-D L +M) J (»'n%^j j- (13) 

The cases denoted by a,b,c,d,e correspond to the cases 2, 4, 11, 14, 16 in the notation of |], ^|. 



Case d, L odd: The ground-state is twofold degenerate in this case. 

<^> ± 4((-l)'cot«^ ± l). (U) 



(<Q = ^-[***^) ■ (15) 



Case d, L even 

-iy f . (2j-l )7rV 1 

• Case e, L odd, 5 = 0: The same as given for case d, L even. 

• Case e, L even, 5 = 0: The same as given for case d, L odd. The ground-state is 
twofold degenerate in this case, too. 

2.2. Exact results in the continuum limit 

The continuum limit defined by L ^> l,z = j/L fixed of the expressions (|TT])-(|H)|) 
is readily obtained. Furthermore we have been able to compute the profile in the 
continuum limit for case e for general values of S (see [Appendix B|) . 



• Case a and case b 

W> = (-1) J (Lsin(^))- 1 cos 2 (^) +0(1/L 2 ). (16) 

• Case c 

K> = (-l) J (^sin(^))- 1 cos(^) +(-l) L /2L + 0(l/L 2 ). (17) 

• Case d: The same expression as given for case e with 5 = and L even and L odd 
interchanged. 

• Case e, L odd 

(of) = (-l) j cos(5(Lsin(7r2))" 1 + 0(l/L 2 ). (18) 

• Case e, L even 

{o*)± = (-iy cos(5 =F ttz) (L sin^r 1 ± l/L + 0(1/L 2 ). (19) 
2.3. Exact results in the limit L — > oo, j ^> 1 

We have obtained exact results for the profiles of a| in the limit L — > oo, j ^> 1 for 
the five cases given in table 0. For the cases a,b,c,d the result may be read of from the 
expressions obtained on the finite lattice. This can not be done for case e with arbitrary 
values of 5, since the respective result on the finite lattice is restricted to 5 = 0. However, 
in the limit L — > oo, j ^> 1 the calculation for case e can be done for arbitrary values of 
S (see [Appendix li| ). Using a z instead of the parameter 5 (see table []] for the relation 



between 5 and aS) we have obtained 



K) = (-i) J /1 I , 2 , r 1 + -. (20) 
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Table 2. Boundary parameters for the cases where we found exact expressions for the 
ground-state profiles of crj in the limit L — > oo,j 3> 1. The first case corresponds to 
Dirichlet-Neumann (DN) boundary conditions, whereas the second case corresponds to 
Neumann-Neumann (NN) boundary conditions on the free boson. The free parameters 
X and ip are real constants. 



case a+ 


a.- 


/3+ 


/3- 




f 1 


1 








DN 


g & 


e -i(P 




e -i(x+v) 


NN 



The profiles for the cases a,b,c,d are also given by this expression if one replaces a z by 
The result (p0|) has to be compared to (|T0|) , which yields B as a function of a z , i.e. 



^ = 7 ON ■ (21) 

We have obtained exact results for the profiles of crj in the limit L — > oo and j > 1 for 
the two cases given in table |2[J]]. Our computations for case f (see [Appendix C.2| ) yield 



(of) = (-l) j Aj- 1 + • • • (22) 
where A is given by 

A = e l ^2 7 / 12 C G 3 (23) 
and Cq is Glaishers constant fnfl, which can be approximately given by 



C G = 1.282 427 129 100 62 ••• . (24) 
For case g we have obtained (see |Appcndix (J.l|) 



(a]) = A cos (pj-i + (25) 

Comparing the value of C in (§) and the value of A in ( p3|) using ([24]) shows that our 
results (^2]) and (|25|) for 9? = are in agreement with (§). 

At this point we want to draw the readers attention to an interesting by-product 
we have obtained during our computations for case g. It is of major interest for our 
conjecture of the continuum limit of the cxj profiles in the case of Neumann- Neumann 
boundary conditions presented in the next section. We have obtained that the ground- 
state profile of (jj for arbitrary values of tp and x ( see table 0) is given by 

(of) = cos L + ^f^j) f(j, L) (26) 



where f(j, L) does not depend on <p, x and may be computed numerically (see |Appendix 
C.I ). The value of m depends on the value of x an d on the lattice length L. For odd L 
we have m = for — it < % < 7t. If x = 71 the ground-state of H is twofold degenerate. 
The value of m is either m = or m = - 1. For even L the value of m is | for > x > — 71 
and — \ for < x < 71 ■ Here the ground-state is twofold degenerate for x — the value 
of m being either | or — |. 

|| In the notation of |l| |) these cases correspond to case 2 respectively case 9. 
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3. Conjectures 

Our results given so far allow to conjecture the ground-state profiles in the continuum 
limit for more general values of the boundary parameters. We have checked our results by 
computing the profiles on the finite chain numerically for lattice lengths up to L = 700. 
Thereafter we have used the finite size data to extrapolate the profiles for L — > oo and 
fixed values of z = j/L. We have performed these extrapolations for 150 different values 
of the boundary parameters each with 19 (L even) respectively 20 (L odd) values of z. 
The accuracy of our extrapolations varies with the choice of the boundary parameters. 
In most of the cases the relative deviation from our conjectures presented in the following 
has been of the order of 10~ 12 — 10~ 7 (see [Appendix D| for details). 



We have considered the following types of boundary terms: 
o Dirichlet-Dirichlet: 

«± = P± = a z = -1= tan Pz = -j= tan ( 27 ) 

Note that for u = | we recover the boundary terms of case e in table |I| 
o Dirichlet-Neumann: Here we have considered only a special type of boundary terms, 
i.e. 

P + = (3_ ^ (3 Z = a± = a z free (28) 
o Neumann- Neumann: We have studied two types of boundary terms, i.e. 

a z = p z = a± = R a e ±iv (3± = R p e M{x+v) (29) 

respectively 

a z = -f3 z a+ = a-= (3+ = /3_ ^ 0. (30) 

Note that the boundary terms ( p9| ) are identical to the boundaries of case g in table [| 
if one sets R a — Rp — 1, whereas no case in table [1] and table corresponds to 
boundary terms which satisfy (p0|) . In this case our result is based on purely numerical 
computations. Observe also, that this is the only type of boundary terms with diagonal 
and non-diagonal terms at one end of the chain, which is studied in this paper. 

3.1. Dirichlet-Dirichlet boundary conditions 

For this type of boundary conditions the ground-state profile of crj vanishes if the ground- 
state is non-degenerate, since the Hamiltonian commutes with of erf • • • er£. Therefore 
we have restricted ourselves to the ground-state profiles of crj. Our results for the cases 
c,d,e from table [1] and our numerical investigations suggest the aj-profiles for general 
values of a z and f3 z , i.e. 

^ j sm(u + 6-2(u + Qir)z) 2(u + Qn) 

This equation is valid for L even and L odd. The value of Q depends on the values of 
uj and the lattice length L. This dependence is given in table |[ The discontinuities 
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Table 3. Dependence of Q on w and L. 





-7T < U) < -f 


-f < w < 


< < f 


| < UJ < 7T 


L even 


Q = i 


g = o 


Q = 


Q = -i 


L odd 


Q = i 


Q = \ 


Q = -| 


Q = -\ 



are due to level crossings in the spectrum at the respective values of u. At these points 
(a? = for L odd and uj = ±ir/2 for L even) the ground-state is twofold degenerate. 

3.2. Dirichlet-Neumann boundary conditions 

Our numerical computations suggest the following profile of the crj-operator: 

i 

(o-z) = ^^A{-l) L+x ^ [sin ( — ) 5 (Lsm(irz))-* (32) 



V 2 

where A is given by (f23|). 

For two cases corresponding to this type of boundaries we have already considered 
the exact <r|-profiles, namely for case a and b from table |l[ The given expression (|16|) 
suggests universal behaviour with respect to the variation of the value of j3±. Our 
numerical calculations for various values of (3± confirmed our guess. Furthermore they 
yield the dependence of the profile on the second parameter a z . We have obtained 



K> = (- 1 ) j t^ cos2 Iy) ( Lsin (^))"- ( 33 ) 

Note that the form of the profile does not depend on the value of a z . However, observe 
also that the amplitude does. For a z = the profile vanishes, since then the Hamiltonian 
H commutes with crfa^ ■ ■ -a x L (the ground-state of H is non-degenerate ||T|). 

3.3. Neumann- Neumann boundary conditions 

Our exact result for case g from table |2|, in particular equation (0), in combination with 
numerical computations yields the crj-profile for the first type of boundaries in question 
(see ©): 

(<7*) = n 1/4 A(-l) j cos (<p + ( X + 2mir)z) (L sin(Trz))^ (34) 

where the value of m depends on the values of \ and the lattice length L as already 
discussed at the end of section 2.3. The cxj profiles vanish, as long as the ground-state is 
non-degenerate (see |Appendix B| for details). We did not consider the crj profiles where 



the ground-state is degenerate. 

For the second type of boundary terms (13(1) the cr| profile does not vanish. For odd 
values of L our numerical investigations suggest 

(a z .) = (-lyVv^ cos(ttz)(L sin(7r^))- 2 . (35) 
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We have not been able to find an analytic expression for the profile for even values of L. 
For a + = the prefactor in fl35|) diverges. However, in this case the profile is already 
given by ([31]). For a z = the profile vanishes as in the case of the Dirichlet-Neumann 
boundary condition (p8[), since in this case the Hamiltonian H has a higher symmetry. 
It commutes with cxfcrf • • ■ o~ x L (the ground-state of H is also non-degenerate here |L|). 

Comparing (|3SJ) with (^) yields that the leading term of the cx|-profile is zero as for 
the boundary terms (|29|), since according to the long range behaviour of the two-point 
function the bulk-scaling dimension associated to the <t| operator is 1. Since there 
exists no spinless primary field with bulk-scaling dimension 2 [TJ|] we conclude that the 
profile (^) is determined by a secondary field (the profile of a primary field with spin 
vanishes f5|). 

We have checked the cxj-profile only for odd L. It is not affected by the presence 
of the diagonal terms in (|30"D and is already given by ([34]) using tp = x = m = 0. 



4. Comparison with the predictions of Conformal Field Theory 

In this section we are going to consider the behaviour of the profiles we have conjectured 
near the boundaries, i.e. we consider the limit 2 < 1 (respectively (1 — z) <C 1) , in 
order to check for the validity of the equations @ and (]!]). 



4-1. Dirichlet-Dirichlet boundary conditions 

The behaviour of (|3l|) near the left boundary is readily obtained: 

= ^ySu^ + S) / , _ + + /I 2u/*\ , + . . \ + %J 

J j 71 " V V" 71 J J 

where u' = u + Qtt. The behaviour of the profile near the right boundary may be 
obtained by exchanging 5 by —5 and performing the transformation j \— > L + l—j. This 
can be seen from (|2"7|). Note that (|3T)| ) does not agree with equation (|5|), since in (^) 
there appears no term which is linear in z. This problem is avoided by considering the 
profile of the following linear combination of cr| and cr| +1 for odd values of j: 

(-m.j-co t (^)^) = -A( 1 + Q_^) M2 + ...). (37) 



The value of the Casimir amplitude A a b can be obtained from the results in || [TT 
It is given by A a b = — + fu- Hence ( p?D is in agreement with (||). Note that this 
operator does not yield the correct scaling behaviour near the right boundary, since 
the given linear combination of erj and <j| +1 depends on the value of 5, which enters 
the left boundary term in a different way than the right boundary term (see (P7j)). 
Another linear combination can be chosen to obtain the correct behaviour near the 
right boundary. 
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4-2. Dirichlet- Neumann boundary conditions 

The Casimir amplitude for this type of boundary conditions is independent from the 
explicit values of the boundary parameters. Its value is given by A a b — ig @- The 
profiles (|32l) and (|3~3"D vanish near the left respectively near the right boundary and 
therefore the predictions (|5|) and @ are not valid. Near the right boundary we obtain 
from © 

K+i-j) = (-1^-' (l " ^ + " " ") (38) 
whereas ( |3~3D yields near the left boundary 

Both expressions are in agreement with the predictions (fj|) and @. 



4-3. Neumann- Neumann boundary conditions 



From (|34] ) we obtain the profiles of cxj near the left boundary, i.e. 

(crp = (-l)Mj'~4 C 0S<y9 

\ /7T 2 (Y + 2m7r) 2 \ \ . , 

1 - tam»( X + 2m7r)z + ( ^ " -J z 2 + ■■■)■ (40) 

The Casimir amplitude has been shown in |]] to be A a b = —■^+2tt (m + ^) 2 . Note that 
the profile is in disagreement with the predicted behaviour for a true scaling operator 
given in (|5|). In order to obtain the correct behaviour near the boundary, we have to 
consider an appropriate linear combination of o x - and ctJ. At the left boundary we obtain 

7T 2 (x + 2m7r) 2 ^ 2 



(cos^xj - sin^xj) = (-iyAj-± ( 1 + ( 24 - 2 ) * + ' ' ' ) (41) 

which is in agreement with (|D and (|6]). Note that building the given linear combination 
is equivalent to performing a rotation around the z-axis such that <p = in (|29|). 
However, similar to the case of the cr| profile for the Dirichlet-Dirichlet boundary 
condition, this linear combination does not lead to the correct behaviour near the right 
boundary. To obtain the correct behaviour near the right boundary one has to perform 
a different rotation as long as x ^ 0. 

It is interesting to consider the respective profile of the orthogonal linear 
combination of and <rj, which yields 

(siny?crj + cos ipaf) = -(-l) j Aj'^ ((x + 2mir)z H ) . (42) 

Note that this profile vanishes near the left boundary and therefore the predictions (|5|) 
and @ do not apply. 
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5. Boundary bound states 

As already mentioned in the introduction, we have also studied the profiles of certain 
excited states of H . These states are obtained by the excitation of a massive fermion 
with respect to the ground-state of H\ ong . We refer to these excitations as to boundary 
excitations, since we are going to see, that they give rise to a contribution to the profile 
which decays exponentially into the bulk. 

First, we are going to look for massive excitations for general, hermitian boundary 
terms in section 5.1 by studying the secular equation p(x 2 ) = which determines the 
fermion energies and which has been obtained in Q. In section 5.2 we will study the 
magnetization profiles near the boundaries of the states, where one massive fermion is 
excited, for two special types of boundary terms. 

In the case of diagonal boundary terms the energy gaps defined with respect to the 
ground-state energy are in one to one correspondence to the boundary 1-strings found as 
solutions of the Bethe Ansatz equations. The 1-strings define the energy gaps between 
the energy of the reference state and the eigenstates in the one magnon sector. At the 



free fermion point a boundary 1-string appears if 2a z > 1 or 2/3 2 > 1 | I3| . The energy 



gaps corresponding to the boundary 1-strings in the limit L —>■ oo are given by 

^Hte; + ^) e ^\{^ +v ~ w ) (43) 

where by E a and E@ we denote the energy gaps corresponding to the boundary 1-string 
which appears for 2a 2 , > 1 respectively 2f3 2 > 1. 

5.1. Massive excitations 

The fermion energies are obtained by 

2A = ^(x + l/x) (44) 

where the values of x are given by the roots of a polynomial p(x 2 ) |IJ]. Since the energies 
have to be real numbers, the roots of the polynomial lie either on the unit circle or on 
the real axis. We will focus on the roots which lie on the real axis for L — > oo. Due to 
(0) these roots yield fermions with a finite energy in this limit. We have not been able 
to determine whether there exist massive fermions corresponding to roots lying on the 
unit circle. 

We are going to restrict ourselves to the search for zeros satisfying \x\ < 1, since 
if x is a zero of p(x 2 ) so is 1/x. In the limit L — > oo, the secular equation p(x 2 ) = 
reduces to 

(1 + (1 - 2a 2 z - 2a^a + )x 2 - 2a 2 z x 4 ) (l + (1 - 2(3 2 z - 2(3-(3+)x 2 - 2(3 2 x 4 ) = 0. (45) 
Solving for x 2 yields two solutions with \x 2 \ < 1, i.e. 

2 1 

X — 



4a? 



1 - 2a_a+ - 2a 2 z + ^8a 2 z + (2a 2 z + 2a.a+ - l) 2 (46) 
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and 



x 2 



i ; , . 2(3.(3+ - 2% + y/8ff + JM + 2/9-/3+ - l) 2 ] • (47) 

Note that these solutions depend only on the parameters of one boundary at a time. The 
fermion energies are obtained by taking the positive square root of ( pf) respectively ([47]) 
and using (f0J). Hence we conclude that the maximum number of massive excitations 
obtained by our ansatz is two (one per boundary). This result is interesting since 



it confirms an analogous result obtained by Ameduri et al |16| in their study of the 
boundary sine-Gordon theory at the free fermion point. 

If the value of a z or (3 X equals zero, the expression ( f4~6|) respectively (|4T| ) has to be 
replaced by 

x 2 = l/(2a_a+ - 1) respectively x 2 = 1/ {2/3-/3+ - 1). (48) 

If a_ = or a+ = respectively (3_ = or (3+ = the result simplifies to 

x 2 = l/(2a 2 z ) respectively x 2 = l/(2/? 2 ). (49) 

All solutions of fl45|) are only asymptotic solutions of the original polynomial if they 
satisfy the assumption \x\ < 1. In the case of diagonal boundary terms this condition 
coincides with the condition for the existence of boundary 1-strings in the Bethe Ansatz. 
The energies of the one-magnon excitations corresponding to the boundary 1-strings 
given in (|43|) correspond exactly to the fermion energies ([0]) obtained by the roots in 



5.2. Magnetization profiles near the boundary 

We have studied the magnetization profiles for certain states, which differ from the 
ground-state of H\ on % by a boundary excitation, for two types of boundaries. Note that 
the ground-state of Hy ong is at least twofold degenerate. We have chosen it such that 
the excited state in question lies in the (+, +)-sector. Therefore these states correspond 
also to eigenstates of H. First we have considered diagonal boundary terms 

a z = 1/(2/3,) = e^/V2 a± = (3+ = (50) 

for even and odd values of L. Without loss of generality we have restricted ourselves 
to £ > 0. In this case 2a 2 becomes larger than 1 and according to ( |4TS| ) and (53) there 



exists one boundary excitation with energy 2A = cosh£. We are able to compute the 
<j|-profile analytically in the limit L — > oo, j ^> 1 (see |Appendix id| ), i.e. 

(t I ^ I - —A^r 1 + 2(e 2 ^ - l)e~^ (51) 
J ix cosh 4 

where |£) denotes the state obtained by the excitation of the massive fermion. Observe 
that the correlation length l/(2£) is related to the excitation mass 2A = cosh£. The 
profile near the right boundary is 
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as expected. The U(l) charge, which is given by the projection of the total spin onto 
the z-axis, of the state | £) is for L even and | for L odd. 

Notice that the exponential part in (|51|) is not present if one considers the respective 
ground-state profile (see (|20|)). We conclude, that adding a massive excitation to the 
ground-state yields an additive contribution to the profile which decays exponentially 
into the bulk. 

The second type of boundaries we have considered are non-diagonal and given by 

a± = a/ cosh £ a exp £ a f3± = ^ cosh £g exp ^ a z = f3 z = 0. (53) 

We restricted ourselves to even values of L, since in this case for odd values of L the 
excitation of a fermion with respect to the ground-state of H\ ong does not lead to a state 
which lies in the (+, +)-sector and hence does not correspond to an eigenstate of H (see 
Appendix C.3| for details). Since the cx|-profile vanishes exactly for the boundaries (|53|) , 
we have studied the <7?-profile in this case. In contrast to the cr|-profiles given in (j5~l|) 
and (0) our results for the crj-profiles we are going to present are based on numerical 
computations for L = oo in the special case ^ = — £ a and for L = 800 in the general 
case. 

In [|T] we have shown that the crj-profiles for the boundaries given in (^) may be 
computed in terms of determinants of j x j matrices. In the case £p = — £ Q it is possible 
to compute these matrices analytically in the limit L — > oo (see |Appendix C.3| ). Thus 
the numerical study of the profile near the boundaries is easier than in the general case. 
Here one boundary excitation appears with energy 2 A = cosh£ Q (see fl4~8| ) and (fP|)). 
For £ a > we have obtained by numerical extrapolation from data for j = 1 . . . 100 

(£a | a] I Q « {-\y +l An (l - 2e(- 1 ^e" 2 ^) (54) 

and 

(a i i u « (-i)^r^ (55) 

where A is given by (^) and | t; a ) denotes the state obtained by the excitation of 
the fermion with energy 2A = cosh^Q,. Again the exponential term is absent in the 
respective expressions for the ground-state profile. Note again the relation between the 
correlation length and the excitation mass. 

Inspired by this result we have also studied the case £ a ^ If £, a ,£,f3 > there 
exist two boundary excitations with energies 2A Q = cosh£ a and 2A^ = cosh<^g in the 
limit L — > oo. We have considered the respective magnetization profiles on the finite 
chain numerically (L = 800). Exemplarily the figures 1 and 2 show the profiles we have 
obtained numerically for the state | £ a ) near the left and the right boundary, where we 
have chosen £ a = 0.075, ^ = 0.05. 

Our results suggest that as long as £ Q ^ £p the profile for the state | £ Q ) is still 
given by (|54]) and (|55|) (the deviation from the expression in (^4|) for large j in figure 1 
is expected due to the presence of the right boundary). The profile for the state | £p) 



15 



10 



0.01 



0.001 



0.0001 



numerical data 

2 e 0. 075-0. 15j 
2 c -0- 075-0. 15j 



15 



20 



25 



30 



35 



45 



J 



Figure 1. The exponential part of the crj-profile of the state | £ Q ) near the left 
boundary for the choice £ Q = 0.075, £p = 0.05 in ([33]). The numerical data (+) have 
been obtained for L — 800. The two lines show the behaviour according to (|54|) for 
even respectively odd values of j. The deviation of the numerical data from the straight 
lines for large values of j is an effect of the finite lattice. 



(-1) J '<6* | crf+l-j I £a> 




100 



Figure 2. The cj-profilc of the state | £ Q ) near the right boundary for the choice 
£, a = 0.075, £/3 = 0.05 in (53). The numerical data (+) have been obtained for L = 800. 



The straight line shows the behaviour according to (B5J) 
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(obtained by the excitation of the fermion with energy 2A ( g = cosh^) near the left 
boundary is just 

whereas at the right boundary we have 

(C/3 I ; , | fc) « (-l) i+1 4f* (l - 2e<-^e-^) . (57) 
However, if £ Q = ^ = £ our computations suggest a different behaviour, i.e. 

(Ca I *J I Q = I ^ I &> ~ (-l^-ie(- l ^e-^ (58) 
near the left boundary and 

(£a I vl+i-j I £<*> = (6 I I fc> « (-l^-ie^V^ (59) 

near the right boundary. In contrast to the profiles we have obtained so far, these profiles 
are left-right symmetric. Furthermore the leading term proportional to does not 
appear, unlike in equations fl54|) and ( f>% . 

Note that the Hamiltonian with boundary conditions (|)3|) conserves parity if 
£<* = Our numerical computations yield that the energies 2A a and 2A^ are non- 
degenerate on the finite chain. Only in the limit L —>■ oo the degeneracy appears. 
Hence the parity of the states | £ Q ) and | is well-defined. This explains why the 
profiles reflect the symmetry of the Hamiltonian. 



6. Conclusions 



In this paper we have studied the one point functions of the <r| and crj operator for the 
XX-chain with hermitian boundaries, which in the continuum limit corresponds to a 
free boson field on a cylinder with Dirichlet-Dirichlet, Dirichlet-Neumann or Neumann- 
Neumann boundary conditions ||. 

No results for the magnetization profiles in the ground-state have previously been 
obtained on the lattice. The only known results have been found by means of the 
bosonization technique |12| and were restricted to the limit L — > oo and j ^> 1 (see 
(§) and fllOl)). The expression fllQf) is valid only for small values of the longitudinal 
boundary field, i.e. for small values of a z in our notation. 

We have obtained the exact ground-state profiles of the cxj-operator on the finite 
chain for the boundary terms given in table 1. The respective profiles are given in 
(|TT|)- (P~5|) . The continuum limit of the profiles is given in (16)-(19). 

The aj-profiles for the cases from table p] in the limit L — > oo, j ^> 1 are given in 
(p0|) . Exact results for the aj-profiles in this limit have been found for the cases in table 
2 and are given in (|22|) respectively (p5|) . Our results obtained on the lattice agree to 
the expressions (§) and (|T0D which have been found using the bosonization technique 
|12|| . Note that in contrast to fllOD our result (|20| ) is not restricted to small values of the 
boundary field. 
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Supported by numerical computations, these exact results have led us to conjectures 
for the profiles in the continuum limit for more general values of the boundary 
parameters (see (|2?1)-(|30|)). For the Dirichlet-Dirichlet case we have studied the most 
general boundary terms fl27j). The cr|-profile for these boundary terms is given in (pl|) 
(the <rj profile vanishes exactly, if the ground-state is non-degenerate). 

We also considered a type of boundary terms, which correspond to the Dirichlet- 
Neumann boundary condition (see (128|)). The o x - profile we conjectured is given in (|32|) . 
The profile of the <r| operator is given in (|33|). 

The boundary terms given by (^) and (|30|) correspond to Neumann-Neumann 
boundary conditions on the boson field. The aj-profile for the first type of boundaries 
(f29l) is given in (|34]) . Here the <r|-profile vanishes exactly, if the ground-state is non- 
degenerate. For the second type of boundary terms (|30|) the profile of the cr| operator 
is given in (p5|) . As in the case of the Dirichlet-Neumann boundary condition the crj 
profile is not affected by the presence of the diagonal boundary terms. 

The expressions we have obtained for the profiles in the case of the Dirichlet- 
Dirichlet respectively the Neumann- Neumann boundary conditions do not coincide with 
the predictions (|j) and (|6|) conformal field theory makes for the profiles of scaling 
operators near a boundary. Here one has to consider certain linear combinations of 
spin-operators to obtain the correct behaviour. These combinations and the respective 
behaviour near the boundary are given in (|37| ) and (|40D. 

In this paper we have also discussed the appearance of massive excitations, i.e. of 
fermions with a non-vanishing energy as L goes to infinity. Such a fermion appears 
if the reciprocal value of one of the expressions ( pESD or ( pE7j ) becomes larger than 1. 
These expressions reduce to (|^) respectively ( f£9|) in the case of purely non-diagonal 
respectively diagonal boundary terms. 

We have computed the profiles near the boundaries in the limit L — > oo, j ^> 1 
for certain states with differ from the ground-state of H\ ong by a massive excitation for 
boundary parameters satisfying ([5(]) or ([531) . For the first type of boundary terms there 
exists one massive excitation, whereas for the second type there might appear two of 
them. The profiles near the boundaries are given by (pT|),(p2|) respectively by (|54])- (|57|) . 
If the value of £, a equals (see (f)3|)) the profiles are given by (|5|) and (|59|). 
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Appendix A. Diagonalization of H\ ong and the projection mechanism 

In this appendix we are going to recapitulate the method we used in [Ij] to diagonalize 
H\ ong given in (§) and the projection mechanism, which yields the eigenvalues and 
eigenvectors of the Hamiltonian H in ([!]). This is necessary in order to fix our notation 
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we are going to use during our computations for the profiles of er| and crj which are 
presented in Appendix B respectively in Appendix C. 
In its diagonal form H\ ong reads 

L+1 L+1 

#long = J2 2AkbkClk ~ J2 Ak ( A ' X ) 

k=0 k=0 

where the bk and satisfy the anti-commutation relations of fermionic creation and 
annihilation operators. The factor 2 in equation ( A.l]) is just an remnant of our notation 
introduced in The values of the A& are determined by the eigenvalues of a skew- 
symmetric (2L + 4) x (2L + 4) matrix M, i.e. 

M(0±) = ±A*(*f). (A.2) 

The explicit form of M is given in Jl|. The eigenvalues Ak can be obtained by the zeros 
Xk of a polynomial p{x 2 ), which has also been given in jL], i.e. 

A fc = \{x k + 1/xfc). (A.3) 

We have also shown in |3J how to compute the eigenvectors (<f> k ) as a function of 
the Xk and the boundary parameters a±,/3±,a z and (3 Z . For certain choices of the 
boundary terms the Xk and the (0^) can be computed exactly on the finite chain, 
since the polynomial p(x 2 ) factorizes into cyclotomic polynomials 0. For the boundary 
parameters given in table [I] and table || this is the case 0J. Otherwise the Xk and the 
(0 fe ) may be computed solving ( A.2|) numerically. 

The eigenvectors (0^ ) can be used to express the Clifford operators rj 1 , which are 
defined by 

in the terms of the creation and annihilation operators bk and ak, i.e. 

L+1 
k=0 

where the (0^) satisfy the orthogonality relation (^fc )j = ^ fc - Equation (|A.5| ) 

allows to write the spin operators in terms of the &&,afe and the eigenvectors (0^). 
Therefore the knowledge of these vectors is a crucial point in view to the computation 
of the magnetization profiles. The vacuum state is defined by 

ak\vac) = V7c. (A. 6) 

Note, that the vacuum state itself is not element of one of the four sectors 
(+, +), (+, — ), (— , — ) or (— , +). However, in view to applications to the Hamiltonian 
H we are interested into the eigenstates of H\ ong which lie in the (+, +)-sector. Hence 
we have defined the states 

|^) = i=(l±6o)l^c) (A.7) 

% They correspond to the factorizable cases 2,4,9,11,14 and 16 in the notation of Q. 
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where bo denotes the creation operator of the spurious zero mode which is always 
present in the spectrum of H\ ong . The associated eigenvectors (0q ) are given by 

1 i 1 i 

O + = (0, -, 0, • • ■ , 0, 0) O = (0, -, 0, ■ • • , 0, ~, 0) (A.8) 

where we used the notation = , (0±)+, (0^, . . . , 

The states [v^) satisfy 

= ±\v ± ) «Ux = ±n\v*) (A.9) 

where rf = 1. If 77 = 1 the state |u + ) is element of the (+, +)-sector and so are all 
the states obtained by an even number of fermion excitations, where we disregard the 
spurious zero mode. If 77 = — 1 all the states obtained by the excitation of an odd 
number of fermions with respect to the state \v~) are elements of the (+, +)-sector and 
hence build up the spectrum of H, where we again have to disregard the spurious zero 
mode. In the case of hermitian boundaries the value of rj may be computed via |TJ 

Note that ( |A.10j) is of no use if there exist additional zero modes on top of the spurious 
zero mode. 



Appendix B. Calculation of the cr| profiles 

In this appendix we are going to explain how to compute the <j|-profiles and how we 
have derived the ground-state profiles for the boundary terms in table |l|, i.e. equations 
(PI)- (POP. We will also derive the profiles ([51]) and ( |52| ) associated to the excitation 
of the massive fermion appearing for the boundaries (0). Note that these boundaries 
coincide with the boundary terms of case e in table [I| if one identifies e^ with tan(| + 1). 
Using (|A.5| ), (|A.7|) and (|A.8| ) one can show that the a z - profile for 1 < j < L and 



for any excited state | ipf) = YY n =i bk n I ^ ,± ) with k n ^ is given by 

r 

I °] I = ^ac I a] I vac) + i£ [(^JMJj " (KWkJj] 

n=l 

where 

L+l 



{vac 



a] I vac) = -i{vac \ r+r^ | vac) = -i J^(0 fc )t(0+ ). . (B.2) 



fc=0 



Before turning to the profiles for the boundaries in table 0, consider the case cx z — (3 z — 0. 
According to |jj it is always possible to chose the eigenvectors (0^) such that the RHS 
of QB.1Q and ([B.2Q vanishes for < j < L + 1 (see the equations (2.38) and (2.40) 



in ). In this case the profiles for all the states | ipf) vanish identically. Thus the 
profile for the ground-state of H vanishes also if it is non-degenerate. The partition 
functions obtained in show that for L ^> 1 this is the generic case. Therefore we 
considered only cases with at least one diagonal boundary term. If the ground-state of 
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H is degenerate one has also to consider linear combinations of the states | ipf) which 
correspond to ground-states H. 

For the cases given in table [1] the polynomial p(x 2 ) factorizes into cyclotomic 
polynomials. This enables us to compute the eigenvectors (0^) exactly on the finite 
chain, as explained in [I]]. The explicit computations are cumbersome but straight 
forward and will therefore not be given. 

For all cases in table |l| there appears an additional zero mode on top of the spurious 
mode which does not give a contribution to the RHS of (|B.1|) . Hence the profiles for the 



two states \v + ) and b z \v~), where b z denotes the creation operator corresponding to the 
additional zero mode, are identical. One of these states is element of the (+, +) sector 
and corresponds to the ground-state of H. However we do not have to know which 
one, since we are only interested into the profile (note that ( |A.10|) is not applicable, due 



to the zero mode). Hence in the following we may assume the state \v + ) to be in the 
(+, +)-sector. 

For the b,c and d with L even no complications arise. The ground-state of 

H is non-degenerate and using the expressions for the {4%) in ( |B.2|) leads to a geometric 
series for each case. Performing the computations yields (|11|)-(|13|) and (|T5|). 



For case d with L odd the ground-state is two-fold degenerate, since here two zero 
modes appear on top of the spurious zero mode. The non-vanishing components of the 
respective eigenvectors are given by 

«)o = 1 = ±i (B.3) 

respectively 



-i V 



i for j even 



\[L I 1 for j odd 



where < j < L and {4>~k)f = ^{4>k)j m ^ n this case we computed the profiles for the 
two states \v + ) and b Z2 b Zl \v + ). This leads to ([L4]) where one has to choose the — sign 
for the profile of \v + ) respectively the + sign for the profile of b Z2 b Zl \v + ). 

However, the computations for case e respectively the boundary terms (^Dj) are more 
involved. It is convenient to use the parameter s = e~ 2 ^ where £ has been introduced in 
(pTl|) instead of 8 being used in table [I]. For odd values of L we have obtained 

where 

rv*^ 3 ^ I />■ ^3 *^ / y» ^3 ^ />-» 2 J ~f" 1 I q ( ,y, ^3 I 1 I rjr> ^3 ^" np ^3 ^ rr> ^3 ' \ 

E x k ' x k ^k x k ~r 5 \ x k x k ^k x k ) /r> p\ 

l/s + s-4-l/xl 

with x k = e 1 ^ and N = {L - l)/2. This expression simplifies to a geometric series for 
s — 1. Straight forward computation yields (|T^) as for case d with L even. We have not 
been able to compute this expression for s ^ 1. 

As for case d with L odd, there appear two zero-modes on top of the spurious 
zero mode if we take L even. The non-vanishing components of the corresponding 



21 



eigenvectors for the first zero mode are again given by ( |B.3| ) whereas for the second zero 
mode they read 

1/2 



22 1 3 



-1 



Z2'3 



L(l + s) 

± 



i for j even 
--y/i for j odd 



(B.7) 



±((/>fc)^. Again we computed the profiles for the two 



where < j < L and (4>k)f 
states \v + ) and b Z2 b Zl \v + ). Labelling the profile for the state \v + ) by a — sign and the 
profile for the state b Z2 b Zl \v + ) by a + sign, we have obtained 



ft 



1 - s 
1 - s L 



si' 1 ± 



1 for j even 
s for j odd 



(B.8) 



^/sL l-s L " L(l + s) 

where is again given by ( B.6|) as for L odd, but with N — L/2 — 1. Choosing s — 1 
allows for the computation of the sum in (|B.6|) and yields flTil) as for case d with L odd. 

Although our results on the finite chain are restricted to s = 1, one can show that 
in the continuum limit [z = j/L fixed, L — > 00) qj is given by 

C— iy 1 

(B.9) 



^ 1/s + 1 sin(7rz 
for L odd respectively by 



-IV 



ft 



COtfTTZJ 



(B.10) 



' J 1/s + 1 

for L even. Plugging this into ( |B.8| ) and ( |B.5| ) respectively results in the equations 
(p"8|) and (0). Note that the exponential contributions in ( |B.5| ) and ( |B.8|) vanish in the 
continuum limit. 

The profile in the limit L — > 00, j fixed, leading to may also be computed. 
In this limit the quantity Sj/L appearing in ( |B.5| ) and ( |B.8| ) can be represented as an 
integral, i.e. 



c°° = lim I 



cos((2j - 3)i) - (a + 1) cos((2j - l)t) + s cos((2j + l)t) 
1/s + s - 2cos(2t) 



dt. (B.ll) 



This integral can be calculated exactly in terms of hyper geometric functions. We 
obtained 



2 

7T 



-IV 



2j - 1 



7T 



and 



2 

7T 



-IV 



2j 



7T 



s 

-IV' s 
j 1 + s 

(-l)'(l -a) 

" (-1) J g 

J 1 + s 



F(l,l/2 + j;3/2 + j;-l/s) 



2j + l 

+ o(r 2 ) 

F(1,1/2-j;3/2-j;- S ) 



for s > 1 



for < s < 1 



+ 



Trfl - s) 



J-1/2 



+ 0(j 



-2n 



(B.12) 
(B.13) 

(B.14) 
(B.15) 
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For the asymptotic expansions of the hypergeometric functions see |T? |. Using these 
results in ( |B.8| ) and ( |B.5| ) gives for even and odd values of L 



(-iy 



■j + 



(B.16) 



7T COsll£" 

Using e^ = \/2a z yields (|20|) . Note that the exponential contributions in (|B.8| ) and ( p.5| ) 
(which, in the limit in question, are present only for < s < 1) are cancelled by the 
exponential part in ( B.15 ). 

We have also considered the profile for a state where the massive fermion is excited, 
i.e. | £) = b^b z \ v + ) for L odd respectively | £) = b^b zi \ v + ) for L even (there are 
two zero modes for L even, where the one labelled by z\ gives no contribution to the 
profile ( p3. 1|) ) . Using (|B.l| ),( |B~T6l ) and the explicit expressions of the eigenvectors (</>^) 
corresponding to the massive fermion, i.e. 

. / i _ e-^ » 1/2 



-1 



-2€0'-i) 



(B.17) 



where < j < L + 1 and ((f)j)f = ±{<f>pf, yields (0) for f > (s < 1). 

Note that for ^ < (s > 1) and j fixed the vector components (|B.17| ) vanish as L 
goes to infinity and the profile for the excited state | £) equals the ground-state profile 
given in ( |B.16|) . By symmetry (see (50|)) this yields the profile near the right boundary 
for the case £ > (s < 1), i.e. equation (|52|). 

The £7(1) charge, i.e. the projection of the total spin onto the z-axis, of the state 
| £) can be computed using the explicit expressions of the (0 fe ) and summing over j 
in (p.l| ) and ( p.2| ) before summing over k (in this way the computation can be done 
exactly on the finite chain, independently of the value of s). This yields for L even 
respectively \ for L odd as stated in section 5.2. 



Appendix C. Calculation of the cr|-profiles 

In this appendix we are going to derive our exact results for the crj profiles for the 
two cases from table |2|, i.e. equations (^),(^) and (|26|). Furthermore we are going to 
explain how we have obtained the profiles (|5^ ) and ( |55| ) corresponding to the boundary 
excitation for the boundary terms in (^3j) with ^ = — £ a . 

We have shown in [1|] how to compute the expectation values of the aj-operator 
in terms of pfaffians. Here we will only consider cases where diagonal boundary terms 
are absent (see table |] and (|53|) ). This enables us to represent (crj) in terms of the 
determinant of a j x j matrix Q, i.e. 



O": 



■)=/« 



Oj 



(T 2 + rx + ) (r+rf) (r+r 3 + ) 



(C.l) 
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where fa = — i for odd values of j and foj = 1 for even values of j. The expectation 
values of the basic contractions of pairs are given by 

fc=0 

The profiles for the excited states 



r 

r> = n^> Hr > ( C - 3 ) 



j'=i 



are also given by the RHS of flC.l|) if one exchanges the (r^Tf) in flC.l|) by 

1=1 

We will start with our computations for case g from table [| and will then turn to case f 



from table Afterwards we will consider the boundary terms (53) . Again it is possible 



to solve the eigenvalue equation (|A.2 ) exactly As for the case of the cx|-profiles we 



will not give the explicit expressions for the eigenvectors ((j>t)- They have been obtained 



as explained in [|TJ []. 
Appendix C.l. Case g 

We will first show, how to derive equation (|26|) . Afterwards we will compute the function 
f(j,L) in @ in the limit L -»• oo,j > 1. Using @ this yields @. 

It turns out, that for our purpose it is convenient to consider the Hamiltonian Hi ong 
in a different basis, i.e. we consider the Hamiltonian 

-^long = UnHx^U^ 1 (C.5) 

where U n is given by 

^(j;)4(-'->(m-V) 

and 

r = -p(,^). (c.7) 

The Oj profile for any state | g) in the original basis is related to the magnetization 
profiles for the state | g') — U \ g) via 

y + mr A . / y + nn \ 

+ L + 1 J J\9 \Oj\9) + sm ( V + L + i J J (9 \Vj\g )• (C8) 

The Hamiltonian in the new basis can still be written in terms of free fermions. Of 
course the fermion energies are not affected by the transformation. The L + l non-zero 
fermion energies are given by 

2A fe = sin f^±^ - (C.9) 







1 ^ 1 i 


]) = COS ( 







L+l 2 L+l 
+ The boundary terms ( |53|) correspond to case 10 in the notation of fl 
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where 1 < k < L + 1. Computing the eigenvectors m the new basis yields 

that they are independent from the values of \ and (p and can be given in terms of the 
vectors (0^ )(0, 0) in the original basis, i.e. for n > we have found 

(^)'(X, <p) = (^ + J(0, 0) for 1 < k < L + 1 - n (CIO) 

(p) = (0^ +n _ (L+1) )(O, 0) for L + 1 - n < k < L + 1 (C.ll) 

whereas for n < 

(0f = (0f_| B |)(O, 0) for |n| + 1 < fc < L + 1 (C.12) 



(^)'(X^) = (0Z + i + ,-i„,)(O,O) for 1 < k < \n\ + 1. (C.13) 



'fc-|n|> 

Now define for n fixed the state 

I n) = b L+2 _ n . . . b L+1 | v^ n ) forn > (C.14) 

respectively 

| n) = b 1 ...b H | v { - 1)n ) forn<0. (C.15) 

For n = define | 0) =| v + ). The expectation values in the new basis may be computed 
in the same way than in the original basis by exchanging the (4>f) in ( |C.4| ) and ( |C.2| ) 
by the (0^)'. Using ( |C.10| )-( |C.13| ) yields that the aj-profiles for the states | n) in the 
new basis are given by 

(n' | a] | ri) = (v + | cr* \ v + ). (C.16) 

where the RHS has to be evaluated in the old basis and for the choice x = f = 0. In 
the same way one can show that 

(n' | a) | n') = (v + | a] \ v + ) (C.17) 

where again the RHS has to be computed in the original basis for the choice \ = tp = 0. 
Since for \ — ¥ — the Hamiltonian H\ ong commutes not only with o"o and cr£ +1 , but 
also with aftcrf ■ ■ ■ cr^ +1 , the crj profile in ( P.17|) vanishes. 

However, we have been interested in the ground-state profiles. Computing the value 
of i] using ( |A.10| ) with (|C.9|) yields r\ = (— 1) L+1 and hence the state | n) is element of 
the (+, (— l) L+1+n )-sector (the excitation of a fermion changes the sector from (±, e) to 
(=F, e) [0). Which of states | n) corresponds to the ground-state of H may be determined 
by computing the associated energies using ( |C.9|) . 

For L odd the ground-state of H corresponds to the state | 0) =| v + ) if x 7r - F° r 
X = 7t the ground-state of H is twofold degenerate. We obtain two linear independent 
ground-states by choosing the states | 0) and | —2). 

For L even and x > the ground-state of H corresponds to | —1) whereas for x < 
the ground-state corresponds to | 1). For x — the ground-state is twofold degenerate 
and we obtain two linear independent ground-states choosing the states | ±1). Finally, 
using ( |C.8| ) and ( |C.16| ) with f(j, L) = (v + \ aj \ v + ) (for the choice x — — 0) and 



m = n/2 yields (^6|) and the values of m as explained in the text. 
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In order to compute f(j, L) we have calculated the expectation values of the basic 
contractions of pairs ( p.2| ) using the vectors ((j>f)(Q, 0). For i, j < L + 1 this results in 

(T±7f> = (r -r+) = v / 2^(j) (C.18) 
(r?Tf)=g(i + j)-g(i-j) (rfr/) = g{% + j) + g{i - j) (C.19) 

where the values of z are always even and the values of j are always odd. Furthermore 
we have defined 

— i sin(r7r/2) 



9 ^ L + lsin(r7r/(2L + 2))' 
Using elementary operations on determinants in ( JU.lj ) one can show that 

/0\L) = (-l) j 0y v^detQ,-, 

where the matrix elements of the j x j-matrix Qj are given by 

f 7i \ 3 f Oh -91 

(Qi)« 



(2k -21 + 1)tt 

sin 



(C.20) 



(C.21) 



(C.22) 



2L + 2J V 2L + 2 

The boundary behaviour of the profile (i.e. (|25|)) is obtained taking the limit L — ► oo 
of Q J5 which yields 

/ 1 -1-1 



lim Q, 



L^oo 



\ — — — 

\ 2j-l 2j-3 2j-5 



3-2j 
1 

5-2j 
1 

7-2j 



(C.23) 



/ 



The large j behaviour of the determinant of this matrix is already known [FS . With 
( |C.21|) we obtain for j > 1 

/0\oc) = (-l)''Ar* + --- (C24) 

where A is defined by (^). Using ( p6|) this yields (|25|) . The function f(j,L) may also 
be calculated on the finite chain using ( |C.22|) in (|C.21|) and computing the determinant 
numerically. 



Appendix C. 2. Case f 

For this case there exists an additional zero mode in the spectrum. Hence the ground- 
state of H corresponds to | v + ) or to b z \ v~), but we can not decide which one, 
since ( |A.10| ) can not be applied. However, the eigenvectors (4>f) corresponding to the 
additional zero mode give no contribution to (|C.4|) as long as i, j < L + 1 . Hence the 
profiles for states | v + ) and b z \ v~) are identical for j < L + 1. Defining 



f(r) 



cos 



2L + 2 



vr g(r) 



(C.25) 
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where g(r) is already given in ( |C.20| ), we have obtained for L odd: 
{rfrf) = (r -r+) = y/2g(j) 

^ T D = /(« + j) - /(» ~ j) ( T * T t) = 9^ + j) + 9{i ~ J) 
where i is always even and j is always odd. For L even we obtain 
(rtrf) = (r -r+) = V2f(j) 

(r+rr) = g(i + j) - g (i - j) ( T r T +) = f(i + j) + f(i - j). 

The expressions ( p.26|) - (fC.29|) have the same large L limit as the respective expressions 
for case g in ( |C18[ ) and ( |U.19| ). Thus the behaviour of the ground-state profile near the 
left boundary is also given by the RHS of equation ( p. 24 ) which yields (P^). 



(C.26) 
(C.27) 

(C.28) 
(C.29) 



Appendix C.3. 

Here we consider the cxj-profile for the boundaries in (|53f) with ^ = — £ a . In this case 
there appears a massive fermion with energy 2 A = 2 cosh £ Q (see section 5). Here we 
are interested into the effect of the excitation of this fermion onto the magnetization 
profile. For simplicity we use £ = £ a in the following. 

We have restricted ourselves to even values of L. Since r\ = (— the state 
|£) = | v~) (where by we denote the creation operator associated to the massive 
fermion) lies in the (+, +)-sector of the Hilbert-space of Hi ong and hence corresponds to 
an eigenstate of H. For odd values of L the (+, +)-sector is given by the excitations of 
an even number of fermions with respect to the ground-state \v + ) and hence the state 
|£) does not correspond to an eigenstate of H for L odd. 

A lengthy computation yields the expectation values of the basic contractions of 
pairs ( |C.2| ), i.e. 

( T ? T i) = 9{i + j) -9{i-j) 

2isinh(2£)e- (i+ ^ 



for?) 



■ g(i+j +1 
1 L + l +9{t 

V2a + f.E(j) 



e 2 « + 1 V L + l 



+ g(j) + 



where g(r) has been defined in ( |C20| ) and 

L_ 



n=l 



' 2£ f 1 ~ 

e s cos 


2 \ 


[ \L + 


I nn ) 



+ e 2 ^ cos 



q 



X 



cosh(2^) — cos 



2n 

TTi 



-71 



L + l 

-i 



-nix 



e -2(r+i)e _ i 
2isinh(20e"^ 

e -2(L+l)£ _ I 



2 cos 



(C.30) 
(C.31) 

(C.32) 



L + l 



-nrt 



(C.33) 
(C.34) 



For L 



oo we can rewrite j^r in terms of an integral which can be calculated using 



standard methods. In this limit we may write 



7T 



2i /sin((i + j>/2) sin((z - j>/2) 



i+j 



(C.35) 
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2i 

7T 



sin((z-j)7r/2) , sin((z + j>/2) 



+ 



i+j 



i+j-l 
2 



\m p (2m+l)5 



2sinl.(2Ue- (i+i)€ ( ^ 1 ^ - aretan(r' ; ' 



m=0 



2m + 1 



(C.36) 



for/-) = - 



2y/2< 



a a.1 



e 2 « + l 



sin(W2) 



i-l 

2 



2 sinh(20e- j€ ( ^ 1 ; arctan(e 5 ' 

m=0 



(C.37) 



2m + 1 

We used these expressions for the numerical calculation of the determinant in (|C.1| ) for 
j = 1, 2 . . . 100. We have chosen 48 different values of £ in the range between —4.0 and 
4.0. Numerical extrapolation yields 

(v + | <tJ I v+> ~ (-l)Mj- 1/4 (C.38) 

where the numerical value of A is given by the RHS of (H) with an accuracy of the order 
of 10- 13 . 

As already mentioned above here we are interested into the profile of the state where 
the massive fermion is excited with respect to \v~). The eigenvectors (0^) associated 
to the massive fermion are given by 
, — a 4 



:/p (0 







" coshe' 

— V / 2exp((l — j)i)/p for j even 
for j odd 



for j even 

iv / 2exp((l — j)i)/p for j odd 



0. 



« 7L+1 " coshe 



exp((l - L)0/p 



where 



and 



sinh(2£) 

-iy (4>h)j ■ m the case £ < this mode does not contribute to 



the limit L — > oo, since ( |C43j ) diverges, i.e. 

<£ | aj | £) ~ {-lYAj- 1 '* for £ < 0. 
This changes for £ > 0. We assumed the following behaviour of the profile: 

(£ | a? | ~ - ce-^)r 1/4 - 



(C.39) 
(C.40) 

(C.41) 
(C.42) 

(C.43) 
in 

(C.44) 
(C.45) 
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Using (|U"35T) -(|C"3l) and (|C39|)-( [C^3|) in O we computed the profile for j = 1, 2.. .100 



for 24 different values of £ in the range between and 4. The values of c and r/ have 
been computed using numerical extrapolation by considering the quantity 



In 1+ , i ' -lnc-jr? (C.46) 





°7 












Our results have confirmed our assumption ( |C.45| ) and the numerical estimates from 
our extrapolations suggest 

c = 2e (-iyv v = 2 £ (C.47) 

which yields fl54"|). The largest deviation of the extrapolated values for c and rj from 
( |U.47| ) has been of the order of 10~ n . We omit any tables on this. Computing the 
profile for £ > near the right boundary is equivalent to the computation of the profile 
near the left boundary for £ < (see fl53p) which is given by ( |U.44 ). This yields 



Appendix D. Numerical verification of the conjectures in section 3 



We have verified our conjectures ([H])- ([35|) for the continuum limit of the magnetization 
profiles for the boundary terms (|27|) -(|30"1). In order to do so, we have used numerical 
diagonalization to obtain the vectors (0jjr) (see ( |A.2|) ) numerically and used ( |B.1|) and 
( C.l|) to obtain the profiles of cr| respectively crj on the finite chain. Notice that ( |C.1| ) 
is only valid if a z — (3 Z — and can not be applied to compute the profile for the 
boundary terms ([30]). In this case we made use of the determinant representation (10.14) 
in |l| of the pfaffian which determines the profiles. 

We computed the profiles on finite lattices with L = 20, 60, 100, . . . respectively 
L = 21, 63, 105, . . ., where the maximum value of L has been chosen in the range of 
500 — 700. For each type of boundaries in (|27|)-(|30D we have used 30 — 50 different values 
of the respective free parameters. 

The values obtained on the finite lattice have been extrapolated to L = oo for 
z = j/L fixed using the BST algorithm, where z = ^ for L even respectively z = |j for 
L odd with n6N. 

As already mentioned in the text, the relative deviation of our numerical results 
from the conjectures is typically of the order of 10~ 12 — 10 -7 . In general the accuracy 



of our computations depends on the choice of the parameters in ( P7|) - (|30T) and we have 
searched for regions in the parameter space where the accuracy of our computations 
breaks down. We found that this is the case for small values of a± and f3± in (p8|) 
respectively in (|29|) and (|30|) . As an example figure Dl shows the relative deviation 
of the extrapolated crj-profile for the Dirichlet-Neumann boundary condition ( p8|) with 
a z = from the exact expression ([32]) for (3± = v8, P± = a/8/20, (3± = v^/lOO and 
for even values of L, where we have used finite size data from L = 20 up to L = 700. 
The deviations are always in the range of the accuracy of the extrapolations. Hence we 
avoid any error bars. 
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Figure Dl. The absolute value of the relative deviation of the extrapolated profile 
(cj)extrapoi from our conjecture (|||) for the boundaries (|2|). The extrapolated profiles 
have been obtained from finite size data for lattice lengths L = 20, 60, . . . , 700. 



The behaviour illustrated in figure [U>1| is not unexpected, since the vanishing of 
a non-diagonal boundary term implies a sudden change of the respective boundary 
condition on the free boson from Neumann to Dirichlet. For the boundaries in question 
this change implies that the a x - profile vanishes exactly even on the finite chain (the 
Hamiltonian H commutes with a\ erf ■ ■ • cr£ for f3± = and the ground-state of H is non- 
degenerate). Since the profiles on the finite chain are smooth functions of all parameters 
appearing in the Hamiltonian, the deviation of the finite size data from the profile ([32]) 
we conjectured for the case (3± ^ in the continuum limit (which does not depend on 
the value of (3±) increases as (3± goes to zero. Assuming that our conjecture is correct, 
it is not surprising that the accuracy of our extrapolations decreases simultaneously. 
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